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Abstract 

We study discrete time Markov processes with periodic or open boundary conditions and 
with inhomogeneous rates in the bulk. The Markov matrices are given by the inhomogeneous 
transfer matrices introduced previously to prove the integrability of quantum spin chains. 
We show that these processes have a simple graphical interpretation and correspond to a 
sequential update. We compute their stationary state using a matrix ansatz and express their 
normalization factors as Schur polynomials. A connection between Bethe roots and Lee-Yang 
zeros is also pointed out. 


LAPTh-033/15 
June 2015 

^nicolas.crampe@univ-montp2.fr 
^kirone.mallick@cea.fr 
®eric.ragoucy@lapth.cnrs.fr 
^matthieu. vanicat@lapth.cnrs.fr 


1 





Introduction 


The study of systems far from equilibrium has greatly benefited from exact results obtained for 
interacting stochastic particle systems [g Eg EH EH Ha SS] and amongst these, the asymmetric 
simple exclusion process (ASEP) has become a paradigm. In one-dimension, the ASEP can be 
solved exactly and many physical quantities (such as its phase diagram, N-point correlations, 
density and current fluctuations etc...) can be calculated analytically daiaEsiEgiii]. Hence, 
the ASEP (as well as the closely related Kardar-Parisi-Zhang equation) belongs to the class 
of integrable models: this exceptional property explains its solvability, its intricate algebraic 
structure and its relations to many different fields of mathematics (representation theory, random 
matrices, combinatorics...). A new field, stochastic integrability, at the crossroads of probability 
and algebra, is indeed emerging, in which the ASEP plays a distinguished role m- 

Integrability is a fragile property that is easily lost when one tries to generalize or deform a 
model: a system has to respect stringent mathematical criteria in order to be integrable (basically, 
the Yang-Baxter relation in the bulk and reflection identities at the boundaries [29ll45j i. Our aim, 
in the present work, is to study integrable, one species exclusion processes with inhomogeneous 
and non local transition rates. We shall consider both cases of periodic (section [T|) and open 
boundary conditions (section [2]): using the transfer matrix formalism, we define discrete time 
stochastic processes, by stating explicitly and representing graphically the updating rules for 
the dynamics. Then, we determine the stationary states of these models exactly. When the 
inhomogeneities Zi that enter the transfer matrix are set to 1, the stationary measure coincides 
with the (periodic or open) TASEP one. 

For the open system, we construct a matrix product Ansatz (23 E] using the Zamolodchikov- 
Faddeev m and the Ghoshal-Zamolodchikov |33j relations, thanks to the formalism developed in 
da Eg EH]. Knowing the invariant measure, we compute analytically some physical observables 
(such as the local density) that will be expressed in terms of Schur polynomials. Thus, our work 
interrelates integrability. Matrix Ansatz and combinatorial properties of interacting particle sys¬ 
tems and contributes to extend the class of exactly solvable models in non-equilibrium statistical 
mechanics. We also find an unexpected relation between the Bethe roots (solutions of the Bethe 
equations found in dzi) and the Lee-Yang zeros of the normalization function (see section ESI). 

1 Markov process on the periodic lattice 

This section is devoted to the construction of a discrete time Markov process on the periodic 
lattice (the site L -|- 1 is identified with the site 1). After discussing the general framework in 
section m we construct in section E2] a one-parameter family of commuting operators t{x\z) 
(using the integrability of the model). Using these operators t{x\z), we define a discrete time 
Markov process that can be interpreted graphically (see section E3D. The stochastic updating 
rules corresponding to this process are given explicitly in section E31 Finally, the stationary state 
of the process and some of its properties are studied in section 11.51 

1.1 General framework: discrete time Markov processes 

We consider a one-dimensional lattice of size L. Each site of the lattice can either be empty or 
carry one particle; a configuration of the system can be denoted by a L-tuple C = ... ,tl) 
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where we set Tj = 0 if the site j is empty and Tj = 1 if it is occupied. The configuration 
space is given by 0 • • • ® where the component of the tensor space 

^ V 

L times 

represents the site of the lattice. A configuration of the system thus corresponds to a vector 

and |1) = 

The dynamics is stochastic: the system can evolve from a configuration C' to another con¬ 
figuration C with the probability M{C,C'). The Markov property is satisfied: the probability 
M{C,C') depends only on the two configurations C and C' and not on the previous history of the 
system. Let Pt{C) be the probability for the system to be in the configuration C at time t. Its 
(discrete) time evolution is given by 



|ti ... tl) = |ri)(8) - • - where we have defined the basis vectors |0) = 


Pt+i{C)= '£M{C,C')PtiC') + 

C'^C 




C'^C 


PtiC) = '£M{C,C')PtiC'), 

c 


( 1 . 1 ) 


where the diagonal term M{C,C) = 1 — probability of remaining in con¬ 

figuration C between t and t -I- 1. This equation can be rewritten in vector-form as follows: 
defining 


l^t) = 


/ Pti (0,...,0,0,0) ) \ 
Pt{ (0,...,0,0,1) ) 
Pt{ (0,...,0,1,0) ) 


E Pti in,.. 

Ti,...,riG{0,l} 


,tl) ) In ...Tl) 


(1.2) 


V Pti (1,... ,1,1,1) ) 

the master equation (jl.ll) becomes 


iP^+i) = M \Pt). 


(1.3) 


A fundamental quantity is the stationary probability distribution of the process, i.e., a vector 
15) such that M\S) = |5) (the existence of this vector is guaranteed because the entries of each 
column of M sum to 1). If the process is not reversible, that is if we do not have the detailed 
balance property M(C',C)5(C) ^ M(C,C')5(C'), the task of finding an explicit expression for |5) 
can be very hard. However in some cases, in particular when the Markov matrix is integrable, 
this stationary distribution can be calculated exactly. 


1.2 Inhomogeneous transfer matrix 


We now use the theory of integrable systems to construct a family of operators, tix\z), called 
transfer matrices, that will allow us to define stochastic updating rules. In this section, we discuss 
only about the integrability of the transfer matrix. Its Markov property will be discussed at the 
end of subsection oi 

The building block of the transfer matrix is the ii-matrix: it acts on two sites and depends 
on a parameter x, called the spectral parameter. Here, we will consider the following i?-matrix 


/1 

0 

0 

0 \ 

0 

0 

X 

0 

0 

1 

1 — X 

0 

V 0 

0 

0 

1 / 


(1.4) 
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It acts in the space 0 C^, where we choose the basis |00), |01), |10), |11) which respectively 
correspond to the following lattice configurations: (0,0), (0,1), (1,0), (1,1). The i?-matrix 
satisfies the Yang-Baxter equation: 


Rul-) Rnl-) R 2 J-) = R23(-) RiJ-) Rul-) , (1.5) 

\X2j \X3j \X3j \X3j \X3j \X2 J 

where the subscripts denote on which tensor space the matrix R has a non trivial action: for 
instance R 12 = R^ 1, R 23 = 1 0 etc... The i?-matrix obeys the regularity relation: 


R{1) = P 


( 1 . 6 ) 


where P is the permutation operator (Pv ^ w = w ^v). It also satisfies the unitarity relation: 

(X) i?21 = 1, (1.7) 

where R 2 i{x) = Pi 2 Ri 2 ix)Pi 2 - The i?-matrix (11.411 is usually associated with the continuous time 
TASEP [2^. 

From the R-matrix, we build an inhomogeneous periodic transfer matrix as follows 


t{x\z) = trQ 


Rq^l 



Ro,L-1 




( 1 . 8 ) 


The parameter x is the spectral parameter and the parameters zi,... ,zl are called inhomo¬ 
geneities. The homogeneous case (i.e. Zi = 1) was studied in [34] . 

The main feature of the transfer matrix is that, due to the Yang-Baxter equation (ll.5p . it 
commutes for different values of the spectral parameter. This is the key to integrability [29] : 


[t{x\z),t{y\z)] = 0. 


We now introduce the following operator 
M{x\z) = t{x\z)t{zi\z)~^ 

= tro 


Ro,L ( — ) Ro,L-1 { ^ 
ZlJ \ZL-1 


Ro,i ( - 


R2,i{-)---RL,ii-)- (1-9) 

Zl Zl 


Note that we have normalized t{x\z) using t{zi\z), but a different choice t{zj\z), j = 2,3, ...,L 
leads to a similar rotated matrix M{x\z). Obviously M{x\z) commutes with t{y\z) and has the 
same eigenvectors. We choose below to use the “normalized matrix” M[x\z) instead of t{x\z) 
because it allows one to construct easily local Hamiltonians (see e.g. remark below). 


Remark: The operator M{x\z) is related to the continuous time TASEP. Taking the homoge¬ 
neous limit Zi —>■ 1, we are left with M{x) = t(x)t(l)~^, where t{x) = tro{Ro^L{x)... Ro^i{x)). By 
standard computation, we obtain 


L 


M'{i) = Y,Pk,k+i.Rik+ia), 

k=l 


( 1 . 10 ) 
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with the convention L + 1 = 1. The operator 


- P.R'{1) 


(° 

0 

0 


0 

0 

1 

0 

0 

0 

-1 

0 

V 0 

0 

0 

0 / 


( 1 . 11 ) 


is the local jump operator of the continuous time TASEP. Hence is the Markov matrix 

of the continuous time TASEP. Therefore, all the results obtained below for the discrete time 
process will be valid, when 2 ^ —>■ 1, for the continuous time TASEP. 

1.3 Graphical representation 

The i?-matrix and the transfer matrices defined above have a useful graphical illustration that 
we now explain. 

R-matrix. The action of the R-matrix Rij{x/y) can be represented by the vertices given in 
fig. E Each vertex corresponds to a matrix element of the R-matrix. A dashed incoming line 


(00|R,,(1)100) (10|R,,(1)110) (01|R,,(|)|10) (10|R,, (|)|01) (11|R,, (|)|11) 



Figure 1: Non-vanishing vertices associated with R. 

(according to the arrows direction) denotes the vector |0) (or equivalently an empty site), whereas 
a continuous thick line denotes the vector |1) (or equivalently an occupied site). In a similar way, 
the out-going lines (after the crossing point) represent the states (0| and (1| on which we are 
contracting the matrix R. The missing vertices in fig. E correspond to vanishing matrix elements 
of R. Remark that the number of particles is conserved: in every non vanishing vertex, the 
number of incoming continuous thick lines is equal to the number of out-going continuous thick 
lines. 

With this graphical interpretation, we will be able to compute efficiently a matrix element of 
a product of R matrices acting in different components of the tensor product. 

Transfer matrix. From expression (11.91) we can deduce a graphical representation for M[x\P). 
The starting point is the lattice illustrated in fig. [2] that one has to fill according to the matrix 
element one wants to compute. Instead of explaining it in full generality, we take below a concrete 
example. 
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Figure 2: Graphical representation for M{x\z). Since the lattice is periodic, the two double slash 
should be considered as linked. 

As an example we take L = 4 and use the graphical interpretation to compute the transition 
rate between the initial configuration (1,1,1,0) and the final configuration (1,1,0,1). The initial 
(resp. final) configuration fixes the form of the incoming (resp. out-going) external lines (dashed 
or thick) as in fig. [3l Then, we look for drawings of the form given in fig. [3] where the remaining 



Figure 3: Starting point for the computation of (1101|M(x|2;)|1110). 

thin lines have to be replaced by thick or dashed lines in such a way that the weights (as given in 
fig- [I]) of all the vertices do not vanish. The total weight of a given possible drawing is then the 
product of all these weights. It is easy to see that there are only two possible drawings, given in 
fig. [Utogether with their corresponding weights. Finally the weight of (1101|M(x|z)|1110) is the 
sum of the weights of the possible drawings. 

Using this graphical interpretation, we are able to compute all the possible rates between 
any two configurations. We remark in particular that the number of particles is conserved (as 
mentioned previously each non-vanishing vertex preserves the number of particles). Therefore, 
we restrict ourselves to a given sector with a fixed number of particles. We can also show that all 
the rates starting from a given configuration (with at least one particle) sum to one which proves 
that M{x\z) can be used as a Markov matrix in the master equation (II.3p . One has to impose 
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also 


( 1 . 12 ) 


0 < — < 1 and 0< — <1, i = l,2, ...L, 

Zl Zi 

so that the probabilities are positive and less than 1. The sector with no particle is special; its 
dimension is 1 and the matrix M{x\z) is reduced to the scalar 1 + — xjzi). Therefore, it 

cannot be considered as a Markov matrix in the empty sector. From now, we consider only the 
cases with at least one particle. 



Figure 4: The two different drawings involved in the computation of the transition rate 
(1101|M(x|z)|1110) with their respective weights. 


1.4 Stochastic updating rules 

The Markov process given hy M{x\z) can be interpreted as a discrete time process with sequential 
update. The configuration at the time t + 1 is obtained from the one at time t by the following 
dynamics: 

• Particle update: starting from right to left (i.e. from the site L to the site L — 1 and so 
on), a particle at the site i jumps to the right on the neighboring site with a probability 
\ — Zi!z\ provided this site is empty. The particle does not jump with the probability Zijz\. 
We remind that we are on a periodic lattice, so that the site on the right of the site L is 
the site 1. Note that a particle located on site 1 does not move. 

• Hole update: once the particle update is done, one performs the hole update. Contrarily 
to the particle update, we do not necessarily start and finish at the sites 1 or L. Let r be 
the site number of a particle (we recall that we restrict ourselves to the case with at least 
one particle). Starting from the site r, we go from left to right up to the site r — 1, using 
periodicity and knowing that the site on the right of the site L is the site 1. A hole at the 
site i jumps to the left on the neighboring site with the probability 1 — xjZi-\ provided the 
site is occupied. The hole may stay at site i with probability xjZi-\. By convention, we set 
ZO = ZL- 

As mentioned previously, we would like to emphasize that, due to the inhomogeneities in the 
transfer matrix, the rates depend on the site where the particle or the hole is situated. All the 
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probabilities are positive and less than 1 thanks to (I1.12I1 . Let us also mention that for the 
homogeneous case, the update simplifies. Indeed, the first step becomes trivial: the particles do 
not move. 

These rules are illustrated in hgure[5]for a chain with 4 sites in the configuration (1,1,1,0) at 
time t. We deduce from this figure the different possible rates between the configurations which 




Figure 5: An example of sequential update corresponding to the Markov matrix M(x\z). The 
first line is the configuration at time t and the third line shows the possible configurations at 
time t + 1. The second line corresponds to the intermediate configurations after the update of 
the particles and the hole update still to be done. The label of the arrows provides the rate of 
the corresponding change of configurations. 


correspond to the entries of M(x\z). One gets 

M(x|z)((l, 1,0,1), (1,1,1,0)) = (1101|M(a:lz)|1110) = ("l - —+ ("l - — (1.13) 

\ ^ 3/^1 \ ^ 1/^1 

in accordance with the calculation done in section [L3] (see also figured]). 

Justification of the sequential up-date. The sequential update described above can be 
easily identified when considering the Markov matrix M{x\z) at the special point x = Zj. Indeed, 
we write M{zj\z) = (t{zj\z)^ where C is the cyclic permutation. Then, from the 

explicit expressions 

t{z,\z) = i?,.,_i(^)...i2,.i(^)i?,.^(^)...i?,.,.+,(i) (1.14) 

Zj-l Zi ZL Zj+I 

t{zi\z)-^ = i? 2 ,l(-)^ 3 ,l(-)---^L,l(-) (1.15) 

Zi Zi Zi 

it is easy to see that ^~^t{zi\z)~^ corresponds to the particle update, while t{zj\z)^ corresponds 
to the hole update at x = Zj. 

Let us remark that such a simple sequential update is specific to the totally asymmetric 
exclusion process. For the partially asymmetric case, it would be much more involved. 














1.5 Stationary State Properties 

We now investigate some properties of the stationary state of the Markov chain defined hy M{x\z). 
We construct the steady state in a formal manner, then give an elementary expression for the 
stationary measure. Finally, we obtain some simple formulas for steady-state correlations. 

1.5.1 Construction of the steady state 

The building block of the stationary state is the following vector 

v[z) = ( 1 ) ■ (1-16) 

It satisfies the Zamolodchikov-Faddeev relation m 

Rl2{zi/Z2)vi{zi)v2{z2) = Vi{zi)v2{z2) ■ (1.17) 

The stationary state of the process M{x\z) is constructed from 

|<5) = vi{zi)v2{z2) ... vl{zl)- (1.18) 

Recall that the subscripts denote which component of the tensor space the vector v belongs to. 
To be more precise, since the process conserves the number of particles, there are L independent 
sectors (as mentioned previously we do not consider the empty sector), each corresponding to a 
given number m of particles in the system. The stationary state is hence degenerate: there is one 
stationary state for each sector. Moreover, the exact normalization of the stationary state 

depends on the sector we are considering. Then, each stationary state is given by the components 
of |5) corresponding to the sector and correctly normalized. 

The following calculation justifies that |<S) given in (I1.18p is the stationary state of the system. 
Indeed, for i = 1, 2, ..., L we have 

t{Zi\z)\S) = ... R^,l{-)R^,L{-) • • • Ri,i+i{^)vi{zi)... vl(zl) 

Zi-1 Zi ZL Zi+I 

= Vi{zi)v2iz2) . . .Vl{zl) 

= I5). 

The first equality is obtained using the regularity property of the R-matrix (II.6h whereas the 
second equality is obtained using L — 1 times the property (ll.lip . Clearly, t{x\z) is a polynomial 
of degree less or equal to L in x. It is possible to show (using the graphical interpretation for 
instance) that in the sectors where there is at least one particle, the degree of t{x\z) is in fact at 
most equal to L — 1. Hence we can deduce through interpolation arguments that in these sectors, 
the stationary state is given by the corresponding components of |<S). Note that when we take 
the homogeneous limit Zj —> 1 we recover the uniform stationary distribution of the continuous 
time TASEP. 

1.5.2 Calculation of steady-state averages 

From the knowledge of the stationary distribution, we can calculate various physical quantities. 
We shall see that some observables can be expressed as symmetric polynomials in the inhomo¬ 
geneity parameters zi, ... ,zl. 
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Probability weight of a given configuration. The weight of the configuration (ti, ... ,tl) 
is readily obtained using (I1.18|) : 

L 

S{ti, ... ,rL) = ]^ + Zido^n) • (1-19) 

i=l 


Normalization factor. In the sector with m particles, the normalization factor of the station¬ 
ary state is obtained by summing the weights of all the configurations with m particles 


Zim)^ Y[zi = eL-m{zi,...,ZL), ( 1 . 20 ) 

/C{1,...,L} i^l 
\I\=L—m 

where ei-m is the elementary symmetric homogeneous polynomial of degree L — m. The nor¬ 
malization factor can be written as a Schur polynomial: 

= s^L-m{zi,... where (1.21) 

L—m m 


We remind the definition of the Schur polynomial (of type A) associated with a partition A = 
(Ai,..., Xl) with Ai > • • • > Al > 0: 


sf(zi, ...,zl) 


det 




det ( z- 


L—i 




( 1 . 22 ) 


This expression of the partition function in terms of the Schur polynomial allows us to relate the 
value of in the homogeneous limit (z* 1) with the dimension of the representation vr"^(A) 

of sl{L) labeled by the Young tableau [A]. Indeed, using the Weyl character formula (see e.g. 
[50l [32] for a review), we obtain 








^ - J 


In particular, we have 


Z{m) 


= dim(7r^(l^-”^)) = 


Z1=--- = ZL = 1 

in accordance with a direct computation starting from (ll.20p . 


L 

m 


(1.23) 


(1.24) 


Density. In the sector with m particles, the particle density at site i is obtained by summing 
the weights of all the configurations with m particles, one of them being at site i 

• • • ? • • • ? ^l) 


(1.25) 


{rii) 


L 

We can show that = m as expected. 

i=l 
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Higher correlation functions. The higher correlation functions take also a very simple form. 
The correlations between the sites < ^2 < • • • < is given by 

—m • • • ; — ■ 2 ' 2 i + 1 ; ■ ■ ■ ? — 1 1 -^ 22 + 1 1 ‘ ‘ ‘ 1 — 1 1 '^ 2 ^+ 1 ; ■ ■ ■ ; 

eL-m{zi, 

For i > m, the correlation functions vanish as expected since the number of particles is m and 
the correlation functions for more than m particles has no meaning. 


(1.26) 


• • • '^ig) 


2 Markov process on the open lattice with boundaries 

We now consider a model on a finite lattice with two boundaries that play the role of reservoirs 
from which particles can enter the system or exit from it. In order to construct an integrable 
system, one needs to solve simultaneously the Yang-Baxter equation in the bulk, leading to the 
i?-matrix, and Sklyanin’s reflection equations at the boundaries [l5], which define two boundary 
operators K and K. These operators allow us to construct an inhomogeneous transfer matrix in 
section EH that we illustrate graphically in section 12.21 and interpret as a stochastic process in 
section E31 The stationary state of this process is computed using the matrix product technique 
[271 [8] that allows us to calculate stationary observables (section EH)- Using the explicit expression 
of the normalisation function, we relate Bethe roots and Lee-Yang zeros in section 12.51 


2.1 Inhomogeneous transfer matrix for the open system 

The building blocks of the transfer matrices for the open system are the R matrix defined in 
section oi and the following boundary matrices: 


K{x) 


{a-\-x)x 

xa+1 

1-x^ 

xa+1 



and K{x) 


xb-\-l 

0 


fcfo+l A 

xb I * 

xb-\-l / 


( 2 . 1 ) 


These matrices act in the space and we choose the basis |0), |1) corresponding to the lattice 
configurations (0), (1) respectively. These iL-matrices are associated with the continuous time 
TASEP with boundaries as will be explained below. 

The matrix K obeys the reflection equation m 


Ri2{ -) Ki{xi) R2 i{xiX2) K2{x 2) = K2{x2) Rl2{xiX2) Ki{xi) i?2l(-) 

X2 X2 


( 2 . 2 ) 


The matrix K is related to an auxiliary K matrix, defined as K{x) = tr^ ^A’o(y)-Roi(^ 
which obeys the reflection equation (12.211 associated with the matrix Ruix) = R 2 i{l/x): 



Ri2{ -) Ki{xi) R2 i{xiX2) K2{x2) 


K2{x 2) Rl2{xiX2) Ki{xi) R2 i{ -) . 


The explicit form of K is 


K{x) 


- 1 ) \ 

x{x + h) 
xb + 1 

x{x + 6) / 


(2.3) 


(2.4) 
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We can now construct the inhomogeneous open transfer matrix as follows: 

t{x\z) = tro {Kq{x) Ro,Li-^) ■ ■ ■ -Ro,i (^) Ko{x) Rifl{xzi)... Rl,o{^zl)^ . (2.5) 

It is shown in [22] that these transfer matrices commute for different values of the spectral 
parameter: [t{x\z),t{y\z)] = 0. Note that the derivation of this property cannot be obtained 
directly from the formalism developed in [13], because the matrix R does not satisfy the crossing 
unitarity. We use the operator t{x\z) to define the following discrete time Markov process 

\Pt+i)=t{x\z)\Pt) . (2.6) 

The parameters must satisfy the following constraints 

0 < xzi <1, 0 < — <1 and ax, bx > 0 (2-7) 

to ensure that the entries of t(x\z) are probabilities. We can also show that the entries on each 
column of t(x\z) sum to one, which guarantees the conservation of the probability \Pt). 


Continuous time limit. In the homogeneous limit Zi —>■ 1, we are left with the operator 
t(x) = t(x\z)\^.^.^^. A straightforward computation gives 


- it'(l) = ^ . (2.8) 

k=l 

where the prime (.)' stands for the derivative with respect to the spectral parameter. Then, 
introducing the parameters a = l/(a + 1) and /3 = 1/(6 + 1), we obtain 



-a 0 \ 
a 0 y ’ 


-P.R'il) 



0 

0 

0 \ 

0 

0 

1 

0 

0 

0 

-1 

0 

V 0 

0 

0 

0 y 



0 /3 \ 

0 -P J ■ 


(2.9) 


Hence —1/2T(1) is nothing but the continuous time Markov matrix of the open TASEP. There¬ 
fore, the stationary state of the continuous time TASEP computed in [27| can be obtained by 
taking the homogeneous limit of the one computed here. 


2.2 Graphical representation of the transition rates 

In fig. [6l each matrix element of the AT-matrix is drawn as a left-reflection diagram. The incoming 
lines (according to the arrows direction) before the reflection point denote the vector the matrix 
K is acting on. A dashed line corresponds to the vector |0) (or equivalently to an empty site), 
whereas a continuous thick line corresponds to |1) (or equivalently to an occupied site). In a 
similar way, the out-going lines (after the reflection point) represent the state of the vector with 
which we are contracting to the left the matrix K. 

Similarly, the graphical representation for the matrix K is given in fig. jT] 

We can now draw diagrams to represent the action of the full transfer matrix. The general 
picture is displayed in fig. [8l 
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(o|K(x)|o) mixm m{x)\i) 



{a-\-x)x 

xa-\-l 


1-x^ 

xa-\-l 


0 


Figure 6: Graphical representation of the iC-matrix. 


(0|iF(x)|0) 



1 

xb+l 


m{x)\i) m{x)\o) 



^ xb+l 

Figure 7: Graphical representation of the iF-matrix. 


(l|iF(x)|l) 



1 


{l\K{x)\l) 



xb 

xb-\-l 


As an example we can compute graphically for L = 1 the transition rate (0|t(x|7)|l) be¬ 
tween the initial configuration (1) and the final configuration (0) (see fig. [9|). The sum of both 
contributions drawn in fig. [9] gives 


(0|t(x|z)|l) 


x{a + x){l - xzi) 

(ax + l)(6x + 1) bx + 1 


2 . x{a + zi) 

^ ’ {ax + l){hx + iy 


( 2 . 10 ) 
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Figure 8: Graphical representation of the transfer matrix. 



1 1 

x{a+x){l—xzi) xzi{l — -^) 

(aa:+l)(6a:+l) fej+l 


Figure 9; Graphical computation of the transition rate (0|t(x|2;)|l) for L = 1. The two different 
contributions are represented with their respective weights. 
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2.3 Dynamical rules 

The transfer matrix t{x\z) defines a discrete time Markov process on the finite size lattice with 
open boundaries. The corresponding stochastic dynamics can be described explicitly using a 
sequential update: starting from a given configuration at time t, the configuration at time t + 1 
is obtained by the following stochastic rules 


• Initialisation: 


— The left boundary is replaced by an additional site (the site 0 with inhomogeneity 
parameter 1) occupied by a particle with probability and unoccupied with prob¬ 
ability 


— The right boundary is replaced by an additional site (the site L + 1 with inhomogene¬ 
ity parameter 1) occupied by a particle with probability and unoccupied with 
probability 


• Particle update: starting from right to left (from site L to site 0), a particle at site i can 
jump to the right on the site i + 1 (provided that the site is empty) with probability 1 — xzi 
and stay at the same place with probability xzi. 


• Hole update: once arrived at site 0, we go the other way starting from left to right (from 
site 1 to site L -|- 1): a hole at site i can jump to the left on the site i — 1 (provided that the 
site is occupied) with probability 1 — and stay at the same place with probability . 

• Summation: Then we have to drop the additional sites 0 and L-|-l and to sum the weights 
corresponding to the same final configuration. 

An example of such sequential update is given in fig. [TU] for L = 1. From this figure we can 
compute the transition rates 


(0|t(a;|z)|l) 


(l|f(x|2;)|l) 


ax 


{ax -I- l){bx + 1) 
1 


XZl I 1-) -|- 1 — XZl 


{ax + l){bx -|- 1) 


(1 — XZl)x'^ + XZl ( 1 — 




x{a + zi) 


{ax + l){bx + 1) ’ 


x{x + b) {1 — x^){l — XZl) 
bx + 1 (ax-|-l)(6x-|-1) ’ 



( 2 . 11 ) 

( 2 . 12 ) 

(2.13) 


The equation (|2.12l) is, of course, identical to the expression (|2.10l) derived using the graphical 
representation in fig. [9l 


2.4 Stationary state and observables 

We now determine the stationary state of the inhomogeneous Markov chain, by using a matrix 
product representation, which will be constructed following a method developed in [22]. We shall 
then compute some physical characteristics of the steady-state. 
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Figure 10: An example of sequential update corresponding to the Markov matrix t{x\z). The first 
line is the configuration at time t. The second line represents the possible configurations after 
adding the two supplementary sites corresponding to the boundaries. The third line corresponds 
to the intermediate configurations after the updates of the particles. The fourth line represents 
the possible configurations after the updates of the holes. The label of the arrows provides the rate 
of the corresponding change of configurations. The last line represents the final configurations at 
time t + 1 after the summation step. 
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2.4.1 Matrix Product Ansatz 

The stationary state is expressed as a matrix product. Indeed, the probability of the configuration 
C = (ri,..., tl) can be written as 




1 


L 

{{W\ll{il-Ti)E{zi)+nD{z,))\V)) , 

i=l 


(2.14) 


where the normalization factor is Zl = llW\C{zi)... C{zl)\V)), with C{z) = E{z) + D{z). The 
algebraic elements E{z) and D{z) can be gathered in a vector 

.4(.) = () . (2.15) 

Using this vector, the stationary state can be written in a more compact way: 

|5) = ^((1U|Ai(zi)A2(z2) ... Al{zl)\V)) . (2.16) 

We now use the formalism first introduced in m and developed in [22l [18] to determine the 
algebra generated by the matrix elements of A(z) and by the boundary vectors. The quadratic 
algebra satisfied in the bulk by E and D will be obtained using the Zamolodchikov-Faddeev 
relations; the boundary Ghoshal-Zamolodchikov relations show how the vectors ((1T| and |U)) are 
coupled to E and D. 

The vector A(z) allows us to write easily the exchange relations between E{z) and D{z): 


Ri2{zi/Z 2)Ai{zi)A2{z2) — A 2 ( 2 : 2 ) Ai (zi). (2.17) 

This Zamolodchikov-Faddeev relation ISB defines an associative algebra because the i?-matrix 
obeys the Yang-Baxter equation. Written explicitly it gives four relations 

' [E{zi),E{z2)]=0 
[Diz^),Diz2)]=0 
f Z1 (zi)E(z2) = Diz2)Eizi) 

, D{z,)Eiz2) = ^ [Eiz2)D{zi) - E{zi)D{z2)] 

Let us remark that the third relation is implied by the fourth. 

The relations between the boundary vectors ((1T|, |U)) and the algebraic elements E(z), D{z) 
are given by the Ghoshal-Zamolodchikov relations 


{{W\K{z)A{l/z) = {{W\A{z), K{z)A{l/z)\V)) = A{z)\V)) . (2.19) 


These boundary relations are consistent with the bulk exchange relations thanks to the reflection 
equation. Remark that the boundary vector |U)) (resp. ((1T|) in fact defines a class of repre¬ 
sentations (resp. dual representations) for the ZF algebra, where |U)) (resp. ((VF|) exists and is 
non-vanishing. We prove below that these classes are not empty by providing an explicit (dual) 
representation that belongs to them. 


17 




Written explicitly we obtain four relations 


{{W\E{z) = {{W\^E{l/z) r D{z)\V)) = ^D{l/z)\V)) 

{{W\Eiz) = ^^{{W\{D{l/z) - D{z)) ™ [ D{z)\V)) = !^{E{z) - E{l/z))\V)) 

( 2 . 20 ) 

Let us remark that for each boundary, the first relation is implied by the second one. 

We now prove that the vector |5) is the stationary state of the transfer matrix. Following 
what was done in [22], one picks up the vector Ai{zi) and moves it completely to the right using 
L — i times the Zamolodchikov-Faddeev relation (|2.17|1 . Next, the vector Ai(zi) can be reflected 
against the right boundary vector |y)) using the Ghoshal-Zamolodchikov relation (|2.19p . After 
that, one can push AiiXjz^ to the left through all the Aj{zj) using L — 1 times (|2.17p . Then, one 
can reflect Ai{l/zi) against the left boundary vector ((VF| using (|2.19l) . Finally, one can move 
back Ai{zi) to its initial position to complete the cycle. Thus, we are left with: 

t{zi\z)\S) = \S) . (2.21) 

If we do the previous cycle the other way round we get 

t{l/zi\z)\S) = \S) . (2.22) 

We show below, using the graphical representation, that the numerators of the entries of t{x\z) 
are polynomials of degree less than L + 2. Since the equation t(x| 2 )|iS) = |5) is satisfied for 2L 
different values of a: (x = zi, 1/zi,..., zl, we get through interpolation arguments that for 

L > 3, t{x\z)\S) = |5) for all x. For L < 5 we checked by computer that this equation is satisfied. 

Degree of the numerator of t{x\z): We use the graphical representation given in figures [H 
El and [T] interpreting the last vertex of fig. [1] as in fig. [TTl One can see that any matrix element 



Figure 11: Vertex with non-intersecting thick lines. 

{C'\t{x\z)\C) can be decomposed into continuous thick lines that do not intersect, as illustrated 
in the figure 

The vertices that do not enter a thick line have degree 0 in x. Thus, the total degree of 
the matrix element is the sum of the degree of each thick line, plus the degree coming from the 
boundaries that we look at separately. We first consider lines that are in the ’bulk’ (that do 
not include a boundary). Looking at the weights of the vertices in fig. [H one can see that the 
incoming part of a line always carry a degree 1. Moreover, a line corresponding to a degree n 
must cross at least n exit lines, the first on the left being always of degree 0. This implies that 
the bulk part is at most of degree L. Figure [7] shows that the right boundary does not change this 
counting. The left boundary can add at most a degree 2, as it can be seen on figurelH Altogether, 
this leads to a total degree L -I- 2. 
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Figure 12: Decomposition of a matrix element into thick continuous lines. 


Explicit representation. We now construct an explicit representation of the algebraic ele¬ 
ments E{z) and D{z). We introduce the shift operators e and 5 such that 5e = 1, which are well 
known in the context of the continuous time open TASEP since they enter the construction of 
the original matrix ansatz m- We define 

E{z) = z + e and D{z) = l/z + 5, (2.23) 

and the boundary vectors such that 

((VF|e = a((W| and 5\V)) = h\V)). (2.24) 

It is readily verified that the Zamolodchikov-Faddeev relation (I2.17P and the Ghoshal-Zamolo- 
dchikov relations (j2.19l) are satisfied. From the results of [27j giving explicit forms for e, 6, ((VF| 
and |F)), we deduce that ((VF| and |F)) can be chosen such that ((VF|F)) / 0 which guarantees 
the non-vanishing of |5). 

2.4.2 Steady-state probabilities 

The matrix Ansatz allows us to calculate the stationary probability of any given configuration. 
Let ^ < ji <■■■< jr < L he integers and C(ji,... ,jr) be the configuration (ri,... ,tl) with 
Tj = 1 if z = and = 0 otherwise. We define the (non-normalized) weight of the word with 
^(j'l), • • •, D{jr) at positions ji,..., jV as 

WLijl,...,jr) = ZLXS{C{ji,...,jr)) (2.25) 

= ((W| ... D{z ,,)... D{zj ,)... D{zj ^)... IF)), (2.26) 

where the dots stand for E{zi) operators. 

The same quantity can be calculated using the explicit representation ()2.23l) : 

H^L(ii,...,jr-) = ((lF|...(,5 + ^)...(5 + ^)...(j + ^)...|F)). (2.27) 

The weight WL{ji,... ,jr) computed from the relations (I2.18h and (I2.2nh is proportional to the 
weight WL{ji,..., jr) computed using the explicit representation (thanks to the unicity of the 
steady-state guaranteed by the Perron-Frobenius theorem). Thus, we have 

(2.28) 
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The multiplicative coefficient is obtained by comparing the weights of the empty configuration: 


f{zi,...,ZL) 


{{W\E{z^)...E{zl)\V)) . 

{a + zi)...{a + zl){{W\V)) ’ 


(2.29) 


f{zi,... ,zl) is symmetric under the permutation of the Zi thanks to (12.181) but also under the 
transformation Zi 1/zi thanks to (I2.20p . The group generated by these transformations is 
denoted by BCl in reference to the Weyl group of the root system of the Lie algebra sp{2L). 
The expression of ... ,jr) is given by 


WL{jl,---,jr) = 


ZI...ZL 

Zj, . . . Zj^ 


E'-' E 

P —0 Qr—p — jr—p 


1 


^qr-p 


Qr—p 1 

E 


1 

^Qr—p—l 


92-1 

-E 


11= ji 


1 

Zqi 



By convention when p = r in the first sum, there is no summation over the qi and the formula 
reduces to 6’’ • We also set jr+i = L + 1 in the last product when p = 0. The 

proof of this formula is obtained by induction on the size L, using the identity 


D{zj^)E{zj^+i) 





— i^jr + l + e) + Zj^ + l 


—E{zj^+i) + Zj^+iD{zj^+i). 
Zj^ 



(2.31) 

(2.32) 


Another important step is to compute the normalization factor of the probability distribution. 
We define C{z) = E{z) + D{z) and C{z) = E{z) + D{z). The normalization factors are thus 
given by Zl{zi, ...,zl) = {{W\C{zi).. . C{zl)\V)), and Zl{zi, ...,zl) = {{W\C{zi) ... C{zl)\V)) 
for the explicit representation. Thanks to the property ()2.28l) . we have 

Zl{zi,. . . ,Zl) = f{zi,...,ZL)ZL{zi,...,ZL). (2.33) 

Zl and Zl are symmetric under BCl. 

Our goal is to get an analytic expression for the normalization factor. For this purpose, for 
any sequence of complex numbers u = (ui,U 2 ,...), we define the shifted product by 


(2|U)* = 


[z - Ui){z - U 2 )'--[z - Uk), k>0, 


^1 if A: = 0, 
We have the result (see proof in appendix lAl): 


(2.34) 


Zl{zi, ...,zl) = 


det {{zj\v)^+'^~^ 


- (l/zj|v)^+2-*) 




det - (l/zj)^+i-*) 


with V = (—a, —6, 0,..., 0). (2.35) 
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When a,b = 0, we recognize in the L.H.S. of equation (12.351) . the expression of the Schur polyno¬ 
mial of type C associated to the partition 1^ = (1,...., 1). We remind that the Schur polynomial 
of type C associated with the partition A = (Ai, ..., Al) with Ai > ... A^ > 0 is defined by: 




= s^{zi,...,zl) = 


a=b=0 


det ((.< 


aZz+I— i-\-Xi 






det(2f+^-*-(l/zj)^+i- 

■ 


(2.36) 


As in the periodic case where the normalization factor is linked to the Schur polynomial of type A 
and to the representation of the Lie algebra sl{L) (see (|1.23p and (|1.24p L the normalization factor 
(j2.36p is given in terms of the Schur polynomial of type C and is associated to representation of 
the Lie algebra sp{2L). These observations will allow us to use results of the Lie algebra theory 
to take the homogeneous limit Zj —>■ 1 (see the remark below). 

When a, b are arbitrary, we need to use some generalizations of the Schur polynomial, called 
shifted (or factorial) Schur polynomials. A shifted Schur polynomial of type A, is defined, for any 
sequence of complex numbers u = {ui,U 2 , ■■■) [36], as follows: 


sf{zi,...,ZL\u) = 


det((zj|u)^-*+^*).. 


det (z 




(2.37) 




Similarly, we define the shifted (or factorial) Schur polynomial of type C a^ 


Sx{zi,...,zl\u) = 


det ((zj|u 




- (lAi|u) 


L-\-l—i-\-Xi 




det (^z}+^-^ - (l/z,)^+i-*) 


(2.38) 




Thus, the normalization factor (12.351) is the shifted Schur polynomial of type C associated to the 
partition and to the sequence v = (—a, —b, 0,0,...). For this particular partition, the shifted 
Schur polynomial of type C can be expanded on the usual type C Schur polynomials as (see proof 
in appendix IBI) 


L+l 

Zl{zi, ... ,zl) = S^L^Zi, ..., zl\v) = ^ 

n=l 


_ in 

a — b 


n(.2:r,..., zj^) . 


(2.39) 


Using the explicit representation, we can also determine the particle density at site i from the 
stationary measure: 


ni{zi,... ,zl) = {b + 


1 A ZL-lizi, Zi-i, Zi+i, ...,zl) 


ZiJ Zl{zi,. .. ,zl) 


+ 


(1 -ab) X 


^ ^ ^ Zi—l-\-k{zi,....,Zi—i,Zi-^-i,...,Zi^^;^ Z]^_i_l_k{Zi^2+k: ■ ■ ■ 1 Zl') 


b=0 


a=0 


k=0 


ZlIzi, ...,zl) 


®It is easy to see that Sx{zi, • • •, ZL\vi) is indeed a polynomial in the variables 21 ,..., zl- 
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Remark: In the homogeneous limit Zi —?> 1, • • •, 1) is equal to the dimension of the 

sp{2L) representation associated to the partition 


’li+i 


-41, •••,1) = 


n 


2L + 2 


^ = dim (1^^^-^)) . 


L \ yZy-t-l — Tl 

This expression is obtained through the general formula (see [32] for instance): 


(2.40) 


dim (7r‘^(A)) = JJ 

1<^<_7<L 

Then from ()2.39p . we get 


Aj — 4 T J — ^ Aj + \j + 2L + 2 — j — i 


J - * 


2L -\-2 — j — i 


n 

2=1 


Xi L \ — i 
L + 1 — i 


V /I 1 9 l{o) — dhib) . ^ ^ n + 1 ( 2L + 2 

Zi(l,l,...,l) =-—- with 5 ^(x) = 2^-^ ( 


a — b 


n=0 


,22+1 


in accordance with the results known for continuous time open TASEP m 

L 


= with /.!(,„) = 53 


P 


i _ i 

a (3 


p=l 


2L — p 


2L — p 
L 


X 


pP+i 


(2.41) 


(2.42) 


The equality between these two expressions is ensured by the identity 

gL{x - 1) = hL{x) - ^ (2.43) 

and recalling that a = ^ — 1 and 6 = ^ — 1. 

In particular, when a = b = 0 (i.e. a = (3 = 1), Zl is equal to the dimension of the 
representation 1^, which is the Catalan number C^+i = 


/ 2L + 2 A 

V E+1 ;• 


2.5 Partition function and Baxter Q-operator 


The models considered here have non-diagonal boundary matrices (see (EH)) and cannot be 
solved by the usual integrability methods. In the last few years, various specific techniques have 
been developed for solving such problems, based on the functional Bethe ansatz [lOl |39l [30], 
the coordinate Bethe ansatz [UlMlISIlllH], the separation of variables [311112], the q-Onsager 
approach [T] and the algebraic Bethe ansatz HE]. Recently, a generalization [121 HI] of the 
TQ relations expresses the eigenvalues of problems in terms of solutions of a new type of Bethe 
equations, called inhomogeneous Bethe equations. 

In this section, we shall find an unexpected relation between the ‘partition function’ Zi and 
the Baxter Q operator appearing in the TQ-relations. 

In |I7] . the Bethe equations corresponding to the eigenvalues and the eigenvectors of the 
transfer matrix (|2.5p have been computed using the modified algebraic Bethe ansatz. In this 
context, the eigenvalue of (12.51) corresponding to the stationary state A(x) is given by 


A(x) = X 


- ^L+l 


b + X -r-r XUk — 1 (x^ — 1) 

1 


bx -t- 


fc=i 


Uk — X [bx -t- I) 


n 


1 


{x-Zj){x-—) 


n 


Uk 


k=l 


Uk-X 


(2.44) 
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where Uk are called Bethe roots and are solutions of the following Bethe equations 




Zp) (Uj Zp 
UjZp 


1 ) 


= [u. 


+'>) n 


fe=l 


Ui — 


Uk 


for j = 1,2,... ,L . 


(2.45) 


Let us stress that here the Bethe roots are the ones corresponding to the steady state. The 
Bethe roots corresponding to other states obey different Bethe equations (see or [23] for the 
homogeneous case). We now introduce the following function: 



(2.46) 


This function (up to a coefficient x^) is linked to the polynomial Q of Baxter |2|: its zeros are 
the Bethe roots. The eigenvalue (I2.44p can be written in terms of the Q function as follows 


x{b + x)Q{l/x) 
{bx + l)Q{x) 


(x^ — 1) 
{bx + l)Q{x) 



- 1) 


XZi 


(2.47) 


For the steady-state eigenvector, we know that A(x) = 1 (see section YTM . Hence Q{x) satisfies 


L 

x{x + b)Q{l/x) — (1 + bx)Q{x) = (x^ “ H 

i=i 

This equation, called TQ relation, allows us to compute explicitly the function Q: for a given L, 
it can be shown that equation (I2.48h has a unique solution of the form (I2.46p . 

For the model studied here, the function Q and the normalization factor Z{zi,... ,zi,) are 
closely related. Namely, we get 

Q{x) = Z{zu...,ZL)\^^_y^ = 


where we have used the explicit form of Z given in (|2.39l) . To prove this result, we remark that 
Q{x) given by ()2.50p has the form (|2.46p and we show that it satisfies the TQ relation (|2.48ll . We 
have 

L+l 

x{x+ b)Q{l/x) - {I+ bx)Q{x) = -x ^ ((-x)” - (-1/x)”') sfi+i_„(zi,..., zl) (2.51) 

n=l 

= (x^ - l)sfz,(zi,... ,zl|x, l/x,0,... ,0) . (2.52) 

Finally, we readily check that ( 2 : 1 ,..., z^jx, 1/x, 0,... , 0) vanishes at the points x = Zj and 
X = l/zj which allows us to conclude that (|2.52p is equal to the L.H.S. of ()2.48p . 


(2.49) 

n=l ' 

-) ^5" sfz,_„_j,(2;i,... ,2 ;l) . (2.50) 

p=0 ^ ^ ^ n=0 


- 1) 


(2.48) 
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Relation (I2.49P means that the Bethe roots (zeros of the function Q) are linked to the zeros 
of the steady-state normalization factor in the complex plane of the transition rate a. These 
zeros appeared previously (for the homogeneous case) in [71 [5] as Lee-Yang zeros and allows 
the generalization of the Lee-Yang theory for the phase transition of non equilibrium system. 
Therefore, relation (I2.49P expresses an unexpected relation between two objects arising from very 
different contexts. 

3 Conclusion 

In this paper, we study discrete time Markov processes with periodic or open boundaries and 
with inhomogeneities. The transfer matrices associated with quantum spin chains are used to 
construct the Markov matrix. We provide a graphical representation of these processes and an 
interpretation in terms of sequential updates. We also study their stationary states using a ma¬ 
trix ansatz developed in |47l [22] and express their normalization in terms of Schur polynomials. 
Finally, in the case with open boundaries, we find a connection between the Bethe roots and the 
Lee-Yang zeros. 

We believe that numerous generalizations of these works are possible and that connections 
with other approaches will be fruitful. For the moment, we studied the models on general grounds, 
without specifying the values of the inhomogeneities. However, it would be interesting to take 
particular choices inspired by physical problems. In these cases, we may be able to go further in 
the computations of the physical quantities. For example, one can choose the inhomogeneities 
randomly, according to a certain law, and study the behavior of the density and of the current. 

To compute the stationary state, we use a generalization of the matrix ansatz but other 
methods have been developed to study such problems. For example, the study of qKZ equations 
[28] [9| seems to be closely related to the approach developed here. A common point of both 
methods is the use of different types of symmetric polynomials (or their deformation). Let us 
also mention the combinatorial interpretation of the matrix ansatz developed in m which may 
be generalized to the case with inhomogeneities. 

The method presented here can be used for other models: for example, in a recent work 
[18] . we classified integrable boundary conditions for the open two species exclusion process 
generalizing the model introduced in [49] (see also unitts]). Thus, discrete time inhomogeneous 
two species exclusion processes can be constructed and studied following the lines presented here. 
The generalization to Al-species also seems possible m- 

Finally, the correspondence between the Bethe roots and the Lee-Yang zeros deserves further 
investigation. We do not know if it is valid only for the TASEP or if it is a general feature. 
Such connection may shed light on the Bethe equations and on the Lee-Yang theory for out-of- 
equilibrium models mM- 
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A Normalization factor and Schur polynomials of type C 

The normalization factor can be expressed in terms of a shifted Schur polynomial of type C 

Z{zx, ...,zl) = ... ,2;l|v) (A.l) 


Proof Let us show be induction on L the previous formula. 
• For L = 1, we have 


Z\{^z\) — ((bF| (z\ -\ -h e + |P)) — z\ - \- a -\- b, 

\ zi J Zl 


in agreement with 


sfi(zi|v) = 


(zi + a)(zi + b) - (^ + a) (^ + 


^1-57 


— Zl -I -h a + 6. 


Let L > 1 such that the property holds for L—1. Consider Zl{zi, ..., zl) and {zi,..., zlI'v) 
as Laurent polynomials in the variable zl- It is easy to see that 


Zl{zi, • • •, Zl) — ((IF| { + -—h e + 5 


ZL-\ -h e + (5 ) |P)) 


Zl^OO 


ZL ((VF| {zi-\ - he + (5 ... ZL-i + 


Zl 


ZL-l 


+ e + S]\V)) 


~ zlZl-i{zi,. .. ,zl-i). 

Zl^OO 


We have also Zl{zi, ..., zl) ~ ZL-iizi, ..., zl-i)/zl- Using the determinant formula for 

ZL^O 

, zl|v), we can develop the determinants in the numerator and denominator with 
respect to the last column to show that (zi,..., z^jv) rsj zls^l-i{zi,. .. ,ZL-i\'v) and 

^ Zl^OO ^ 

s9l(zi, • • •, zlW) ~ s9l-i{zi, ..., zl-iW)/Zl- Using the induction hypothesis we deduce that 


s^l{zi,...,zl\v) ~ Zl{zi,. .. ,zl) and s^l{zi, ... ,zl\v) ~ Zl{zi, ..., zl). 

Zi-^oo ^ ZL—>-0 


Hence Zl{zi, ..., zl) and s^l{zi, ..., z^Iv) are Laurent polynomials in the variable zl of degree 
1 in Zl and IjzL and they have the same coefficients in front of zl and I/zl- So we can write 

Zl{zi,.. .,Zl)- S^l{zi,. . ■,Zl\v) = f{zi, . . ..ZL-l), 

with / independent of zl- But Zl{zi, ..., zl) and s^l(^I) • • • j 2 ^l|v) are symmetric under the 
permutation of the Zi's, so / cannot depend of the ZiS and hence is constant. We compute 
this constant by evaluating the previous expression at the particular point where all the z^’s are 
equal to 1. On one hand the value of s^l( 1, • • • , l|v) can be computed using (I2.39P and (|2.40l) . 
On the other hand the value of ... , 1) is equal to the value of the usual continuous time 

TASEP normalization factor, which matches with the previous value (see (12.431) 1. Hence we have 
^£,(1,..., 1) = s^l( 1) • • • ) l|v), so / = 0. This concludes the proof. 
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B Expansion of shifted Schur polynomials of type C 


We show that the type C shifted Schur polynomial (I2.35P can be expanded on the regular type 
C Schur polynomials as follows 


L+l 

sfi(2;i,...,ZL|v) = ^ 
n=l 


gu _ ya 

a — b 


s 


C 

IL+I- 


n 




(B.l) 


where v = (—a, —b, 0,..., 0). 


Proof. Starting from 


det((z,|v)^+^--(l/z,|v)^+^-)^. 

det 


the general idea is to use the multi-linearity of the determinant in the numerator to expand it 
and recognize non shifted Schur polynomials of type C. First we rewrite 


(Zj|v)^+2-* - (l/Zj-|v)^+2-* 


L+2—2 
^3 


L+ 2 -i + (« + ^) 



1 


^3 


+ab 




(B.2) 


Define the matrix M{zi,..., zl) = {{zj\v)^~^‘^ ^ — {1/Zj\v)^~^‘^ ^)ij' Ri{zi,..., zl) be the 
i-th row of M( 2 :i, ..., zl)- The expansion (jB.2ll . we can rewritten as 


Ri{zi, ...,zl) = Pi-i{zi, ...,ZL) + {a + b)Pi{zi,..., zl) + abPi+i{zi, zl), (B.3) 

with Pk{zi,..., Zl) the row vector {zj'~^^~^ — l/z^~^^~’^)j. Remark that Pl+i{zi, ..., zl) = 0. We 
develop det {M{zi,..., zl)) nsing the expansion ()B.3jl and the multi-linearity of the determinant. 
Since we have only L -|- 1 different row vectors Pi, antisymmetry of the determinant imposes that 
we need to choose L distinct row vectors Pi among Pq, ... ,Pl in order to get a non vanishing 
determinant (hence there is just one Pi missing in each non vanishing determinant). Then, the 
expansion of det(M(zi,..., zl)) reads 


where 


L+l 

det(M(zi,... ,zl)) = y^^CnUn{zi,...,ZL), 

n=l 


/ Po{zi,...,zl) \ 


UL+i-k{zi, ...,zl) = det 


Pk-i{zi,. . . ,zl) 
Pk+i{zi,..., zl) 


with k € {0,.., L} 


\ Pl{zi, ..., zl) / 


(B.4) 
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and Cn are some coefficients to be determined. 

Let us first compute Ci (which corresponds to the term with missing Pi in ()B.4p ). There is 
only one possibility to get the row vector Pq: it comes necessarily from the row Ri. Then Pi 
appears in the rows Pi and P 2 , but since we have already chosen Pq for the row Pi, Pi has to 
be taken in row P 2 . One shows by induction that for all i G {0 ,... L — 1}, Pi has to be taken in 
row Pi+i- It is summarized on the following diagram. 

Pi Po Pi P2 

P2 Pi P2 P3 

Pl-1 Pl-2 

Rl 


(B.5) 

Pl-1 

Pl-1 


Each line in (|B.5I1 corresponds to the expansion (IB.31) of a row vector Pj. The coefficients have 
been omitted, but can be easily recovered from (jB.31) . The bold Pj’s in the diagram are the ones 
chosen in the expansion to obtain the term ui(zi,..., zi). The shaded Pj’s are the ones that do 
not appear in 111 ( 2 : 1 ,..., zi). In the case considered here, all the chosen Pj’s correspond to the 
first term in (|B.3I1 and hence come with a factor 1. We deduce that Ci = 1. 

The computation of C 2 (with Pl_i missing) looks similar. There is again only one choice for 
the Pi’s to get U 2 {zi,..., zi). It is given by 


Pl 

Po Pl P 2 



P 2 

Pl P 2 Pi 



Pl-2 

Pl-3 Pl-2 

f}m 


Pl-1 

Pl-2 

mnt 

Pl 

Rl 


Mf 

Pl 


The chosen Pq, ... ,Pl -2 correspond to the first term in ()B.3jl and hence come with a factor 1, 
but the chosen Pi corresponds to the second term in (IB.3P and thus comes with a factor a + h. 
We deduce that C 2 = a + h. 

To compute Cn (with Pi_|_i_„ missing), we draw the diagram 

Pi Po Pi P2 

P2 Pi P2 P3 


Rl—ti 

Rl+I—ti 
Pl+ 2—n 

Rl+S—ti 


Pl— 1—n Pl—ti 

Pl-ii 


Pl+ 2—n 

Pl+ 2—n Ph+'i—n 

Pl+ 2—n Ph+'i—n Ph+A—n 


(B.7) 
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In opposition to the previous cases, there is some freedom in this diagram. Indeed, although the 
first rows Ri,... ,RL+i-n uniquely fix which Pi’s are chosen (written in bold), the row RiJ^ 2 -n 
corresponds to two choices: 

• we choose PL+ 2 -n (which comes with the factor a + b). Then, the remaining rows just 
correspond to the computation of Cn-i- 

• we choose Pi+s-n (which comes with the factor ab). In this case, we have to choose PL+ 2 -n 
in the row R]^^ 3 _n (which comes with the factor 1). Then, the remaining rows correspond 
to the computation of Cn -2 (except that the order of PL+ 2 -n and P^^^-n is inverted, hence 
a factor (-1) by antisymmetry of the determinant). 

Altogether, we get a recursion relation for the coefficients Cn- 

Cn = {a + b) Cn-i — ab Cn- 2 - (B.8) 

It is straightforward to check that Cn = verifies (jB.Sp . Ci = 1 and C 2 = a + b. We 

complete the proof of (jB.ip by remarking that Un{zi,... ,zl )/det — (l/zj)^^^”*) = 

s^L+i-n{zi,. - - ,zl). 
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